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ABSTRACT 

We discuss the gravitational wave background (GWB) from a cosmological population 
of gamma-ray bursts (GRBs). Among various emission mechanisms for the gravita- 
tional waves (GWs), we pay a particular attention to the vast anisotropic neutrino 
emissions from the accretion disk around the black hole formed after the so-called failed 
supernova explosions. The produced GWs by such mechanism are known as burst with 
memory, which could dominate over the low-frequency regime below ~ 10Hz. To esti- 
mate their amplitudes, we derive general analytic formulae for gravitational waveform 
from the axisymmctric jets. Based on the formulae, we first quantify the the spectrum 
of GWs from a single GRB. Then, summing up its cosmological population, we find 
that the resultant value of the density parameter becomes roughly Slew ~ 10~ 20 over 
the wide-band of the low-frequency region, / ~ 10~ 4 — 10 1 Hz. The amplitude of GWB 
is sufficiently smaller than the primordial GWBs originated from an inflationary epoch 
and far below the detection limit. 
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The observation of gravitational waves (GWs) is one of 
the most important missions to explore the trackless parts 
of the Universe. Several ground-based laser interferometers 
(TAMA300, LIGO, and GEO600) are now operating and 
taking continuously the data at a frequency range, 10Hz - 
10kHz, where rapidly-moving stellar objects accompanied 
with strong gravity, such as formations of neutron stars 
or black h oles, are the most promising sources for GWs 



(See iNewl 1120031) for a review). The Laser Interferometer 
Space Antenna, LISA 1 covering 10 -4 — 10~ 2 Hz will be 
launched near future, and more over, the future space mis- 
sions such as BBO 2 or DECIGO dSeto et all200ll) target the 
frequency window around 0.1Hz. For such interferometers, 
the primordial GWs generated during an inflationary epoch 
could be detected. Also, these interferometers are expected 
to be useful in determining the cosmic equation of the state 
llTakahashi fc N akamura 2005]). Thus, low-frequency GWs 
detected via space interferometers may provide a power- 
ful cosmological tool to probe the extremely early Universe 
jMaggiorell200Ct) . One important remark is that several as- 
trophysical foregrounds around deci-Hertz band proposed so 
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far are believed to be weak or resolvable (e.g .,|Ferra ri et alJ 
ll999USchneider et al.ll200ll: iFarmer fc Phinnevll2003D . 

Recently, however, iBuonanno et alJ ll2004) pointed out 
the possibility that a cosmological population of core- 
collapse supernovae could contaminate the inflationary GWs 
at the low-frequency range mentioned above, as a result of 
the GW associated with the neutrino emissions. Their anal- 
ysis suggests that even in a deci-hertz band, one may not 
simply neglect the contributions from the cosmological stel- 
lar objects, especially with the asymmetric emissions of neu- 
trinos. The important remark is that GWs from neutrino- 
driven jets are generated by different emission mechanisms 
from those of the periodic GWs. For an energy flow orig- 
inating from a burst, the GW amplitude suddenly rises 
from zero and settles down into a non-vanishing value af- 
ter the burst. This phenom enon is called burst with memory 
jBraeinskv fc Thornell987l) . It has been argued that a burst 
of neutrinos released anisotropically can generate a burst of 
GW accompanying with the memory effect. A cosmological 
contribution of such GW may therefore become significant 
at the low-frequency band rather than the high-frequency 
regime. 

Among several candidates other than the core-collapse 
supernovae, GRBs are expected to have powerful asymmet- 
ric jets with the enormous emission of high energy neutri- 
nos. It has been already pointed out that relativistic jets 
of matter associated with GRBs can generate strong GWs 
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enough to be det ected by the present laser interferometers 
JSago et all2004h . In fact, a data analysis searching for GWs 
associated with the very bright GR B030329 have been per- 
formed by LIGO group jAbbott et alJl2005h . although they 
reported a null result. 

Another important aspect is that these jets are thought 
to be driven by thermal energy deposition due to the 
neutrino and anti-neutrino annihilation into electrons 
and positrons, which occurs primarily at the polar re- 



gion of the black hole- accret ion disk sy s tem, namely, the 
neutr in o-driven GRBs dPacz ynsk 3 Il990t iMeszaros 



19921: IPopham et al 



3 LL2 




Reed 



iMacFadven fc Wooslevl 119991: lAsano fc Fukuvamal 12000. 
(see, iProea et alj 12003 : iMizuno et alJl2004 iTakiwaki et alJ 



1996 
120051 



iFrver et alJ 



2004; 



se supernova 


3 iBurrows fc Havet 


iMuller et all 


l2004t 


Kotake et al 



, anisotropic energy flows from the GRBs also become 
important sources of G Ws, in addition t o the anisotropic 
jets of the matter itself dSaeo et al.ll2004l) . While a cosmo- 
logical contribution of those GWs may not play a significant 
role due to the relatively rare event rate of the GRBs, a 
quantitative estimate of the shape and the amplitude of 
the spectrum still remains important issue to clarify the 
nature of gravitational wave background arising from the 
memory effect. We then wish to understand how the GW 
memory is generated and how much amount of the GWs 
is released from a cosmological population of GRBs in an 
analytic manner. 

The purpose of this paper is to give an analytic estimate 
of the gravitational wave spectrum from the GW memory ef- 
fect, particularly focusing on the neutrino-driven GRBs. In 
Sec |2| after briefly describing the axisymmetric jet model, 
we derive a general analytic formula for gravitational wave- 
form with memory effect from the asymmetric relativistic 
jets, which is both applicable to the ultra-relativistic and 
the relativistic cases. Using this formula, in Sec|3 we com- 
pute the spectrum of GWs from a single neutrino-driven 
GRB. Then, in SecQ] summing up a cosmological popula- 
tion of GRBs, we estimate the amplitude of gravitational 
wave background (GWB). Uncertainties of the model pa- 
rameters are also taken into account when evaluating the 
GWB, and we find that the resultant GWB amplitude be- 
comes rather small and would be difficult to detect by the 
present and future missions of GW observatories. Finally, 
summary and discussion are devoted to section 



2 GRAVITATIONAL WAVES FROM 
RELATIVISTIC JETS 

In this paper, we specifically treat the GWs from the ultra- 
relativistic jets, particularly focusing on the neutrino-driven 
GRBs as a plausible GW source. In this case, energy sources 
powering the GRBs rely on th e neutrino emission from 
the accreting disk iPiran| [ 2002L for a r eview and refer- 
ences therein). Recently. ISetiawan et all i2004l) performed 
the three-dimensional hydrodynamic simulations and ob- 
tained the time evolution of the luminosity of neutrinos 
emitted from accretion disks around hyper-accre t ing st ellar- 
mass black hole. Using their results. fAlov et all l)2004h per- 
formed the two-dimensional simulations to study the evolu- 



tion of relativistic jets driven by thermal energy deposition 
due to ^F-annihilation. In next section, we adopt their fit- 
ting result for the neutrino luminosity evolution to calculate 
the GW amplitude. Before estimating this, we discuss the 
GW memory and derive a general analytic formula for grav- 
itational waveform from the axisymmetric relativistic jets. 

According to the memory effect of GWs, the GW am- 
plitude h(t) jumps from zero to a non-vanishing value and 
it keeps maintaining the non-vanis hing valu e even after the 
energy source of GWs disappears. lEpstein| lll978l) has de- 
rived a general formula of memory effect produced by the 
radial emission of neutrinos. The detectability of such effect 
was discussed by iBraginskv fc Thornd JlQ87fl through the 
groun d-based laser interferometers. Recently. ISeealis fc Oril 
teOOlI) studied the memory effect generated by a point par- 
ticle whose velocity chang es via gravitation al interactions 
with other objects. Further, ISago et al.l l| 2004h estimated the 
amplitude of GW memory from the GRB jets based on the 
internal shock model. Here, we consider the GWs from the 
neutrino-driven GRB and derive a useful analytic formula 
for the axisymmetric emission of neutrinos. 

Let us introduce the two coordinate systems shown in 
figure^ the source coordinate system (x' , y , z') and the ob- 
server coordinate system (x,y,z). In these coordinate sys- 
tems, the z'-axis is chosen as jet axis, while the z-axis is set 
to a line-of-sight direction. Further, the origins of these two 
coordinate systems are set to the centre of a GRB object. 
Then, the viewing-angle denoted by £ is given by the an- 
gle between z- and z'-axis. For convenience, we assume that 
the z-axis lies on the (x 1 , z')-plane. In this case, the two po- 
larization states of GWs satisfying the transverse-traceless 
conditions become h+ = h xx — — h yy and h x = h xy = h yx in 
the observer coordinates and we obtain h x = (see below). 
Note that the sum of the squared amplitudes, \h+\ 2 + \h x | 2 , 
is invariant under the rotation about the z-axis. The geo- 
metrical setup shown in figure Qyields the following relation 
between the two polar coordinate systems (9, (f>) and (#', <f>'): 

sin 0' cos 4> = sin 9 cos <f> cos £ + cos 9 sin £ , 
sin 9 sin <j> = sin 9 sin (j>, 

cos 9' = — sin 9 cos <f> sin £ + cos#cos£. (1) 

Since the GRBs are cosmological, one can safely ne- 
glect the time lag betwee n the two points near the source 
region, i.e., \t - 1'\ < r/c (iMiiller fc Jankall2004 . Thus the 
amplitude of GWs observed by the distant observer may be 
written as 



, ,_ /s /3 2 sin 2 9 

ifl a cos2 ^ 

1 — p cos 9 



(2) 



where j3 is the velocity of matter in a jet normalized by c, 
and 9 and (j> are related to 9' and 4>' through the relation 
0. The function tp(Qf) represents the depe ndence of the 
amplitude on the direction of the energy flow JSeealis fc Oril 
l200lT) . Note that the counter part of the amplitude, h x , is 
obtained just replacing cos 20 with sin 20 in equation 
which immediately yields h x — by integrating over the 
angle cj>' ■ The quantity dL/dQ! is the neutrino luminosity 
per unit solid angle emitted from the source. For an axisym- 
metric neutrino emission, we consider the two axisymmetric 
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jets whose luminosity is uniformly distributed with respect 
to the opening angle of jets. Denoting the opening angle by 



!, we have 



e 



E 



An 



jet 



T(t) 



e 



(3) 



where the function 0(6*) is Heaviside step function. The 
quantity E is the total isotropic emission energy, which 
is typically E ~ 7?nc 2 ~ 10 51 erg in the case of jets 
from relativistic mat ter, and E v < 10 5 2 erg in the case 
of neutrino emission l|Popham et, alJll999r . In equation ©, 
Afijot denotes the solid angle of two jets given by Af2j e t = 
4n{l — cos(f9 opcn /2)}. The function T(t) represents the time 
dependence of the luminosity of the ejected matter or neu- 
trinos, which will be given explicitly. 

Substituting equation into equation and 

performing the Fourier transformation 2mfh(f) = 
J^° aa h(t)e~ 27V *f t dt, the characteristic amplitude defined by 

hc(f) = fh{f) is given by 

G # 



hc(f) = 



irc 4 r Att 



jet 



-E\W)\ 



(4) 
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Here, the function f (/) is the Fourier transform of the func- 
tion F. The quantity ^ is the angular integral of the function 
ip(Sl'), which depends on the parameters, #o P en,C and /3: 



dQ!ip(Q!) 



2n 



/3 sin 2 £ 



- 2/3X(l + cos 2 £) 



(5) 



+ {/3X - cos £) \J (/3cos£ — A) 2 + (1 — /3 2 )(1 — A 2 ) 
+ (PX + cos £) V(/3cosC + A) 2 + (l -/3 2 )(1- A 2 )| 

with A being cos(f9 opcn /2). In performing the integration, the 
range of the integral denoted by "jet" is restricted to the jet 
region, i.e., < 6' < 9 opcn /2 and -it — 6 opcn /2 < 0' < ir. In 

1) 



the ultra-relativistic limit (j3 

| cos C| 

I *±7I /\ I ; 

* = 



4tt(1 - A) 



sin^ 
A — 



cos 2 £ 



this yields 
for | cos C| > A, 
for | cos C| ^ A. 



(6) 



Figure|2|illustrates the viewing-angle dependence of the 
GW amplitude $/Af2j e t, in the cases of relativistic (/3 = 0.9) 
and ultra-relativistic limit [fi — 1.0). Figure clearly shows 
that the intensity of the GWs generated from an axisymmet- 
ric relativistic jet has an anti-beaming distribution. Com- 
pared to the relativistic case (/3 < 1), the anti-beaming fea- 
ture of ultra-relativistic jet is very sensitive to the opening 
angle, while the viewing-angle dependence becomes fairly 
weak for a typica l value of the opening angle, #o P cn ^ tt/6 
jFrail et al.ll200 ill . This indicates that the GWs from the 
neutrino-driven GRB are uniformly emitted in contrast to 
the axisymmetric emission of the neutrinos itself. 

Analytic expression with JjJ an< i © i s one OI the 
main result in this paper. It describes the gravitational wave- 




jet region 



y.y 



Figure 1. Source coordinate system (x',y',z') and observer co- 
ordinate system (a;, y, z). The observer resides at the distant point 
on the 2-axis. The z'-axis coincides with the symmetric axis of 
the GRB source. Inside the jet represented by dashed lines, the 
luminosity of energy flow is assumed to be uniformly distributed 
with the opening angle 9 open . 
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Figure 2. Viewing-angle dependence of the amplitude. The up- 
per three lines represent the ultra-relativistic cases (/3 = 1.0), and 
lower ones relativistic cases (/3 = 0.9). We take three kinds of the 
opening-angle, 9 opcn = 1/100, ir/12, 7r/6. 



form generated from a (ultra-)relativistic emission of neu- 
trinos and/or matter with opening angle 8 opan and with the 
velocity (3. With the analytic expression, one can further es- 
timate the GWB spectrum in a tractable manner. Before 
doing this, we first apply this formula to estimate the am- 
plitude of GW from a single neutrino-driven GRB in next 
section. 
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3 GRAVITATIONAL WAVES FROM A SINGLE 
GRB 

Provided the time evolution of the luminosity T(t), the an- 
alytic formulae gj with equation (J5J and/or © enable us 
to estimate the frequency dependence of the characteristic 
amplitude, depending on the parameters (£, opcn , E, r). In 
this section, for a simple but a realistic model of the lumi- 
nosity evolution, we adopt the fitting result bv lAlov et alJ 
i2004) to calculate the GW spectrum from a single GRB. 

According to the stat e-of-the-art three-dim ensional hy- 
drodynamic simulations iSetiawan et alj|2004l) . the fitting 
function is given by 



r(f) = 



6tf - tb 



(t/th) 
(t/U 



-3/2 



< t < t b 
tb ^ t < tf , 

tf < t 



(7) 



and T(t) — for t < 0. Here, the time tb roughly corresponds 
to ^ ~ U/3. Note that the three phases of the time evolu- 
tion in equation Q are intimately related to the evolution- 
ary phases of the neutrino emission; i) the mass accretion, 
which makes the luminosity linearly growing, ii) stationary 
phase after stopping the mass fall, and iii) slow decay by 
the termination of neutrino production in the accretion disk. 
These basic behaviors should be fairly insensiti ve to a de- 
tailed modeling of neutrino emission. Although lAlov et alJ 
i2004f) only considered a short- duration GRB in which the 
duration tb and tf were set to 0.01 sec and 0.03 sec, respec- 
tively, we keep to use equation J7J for a long-duration GRB. 
Additionally, note that the observed duration is a few times 
larger than the duration tf due to the tail part of equation 
Q at the late time. The observed duration is typically 0.2 
sec for short-duration GRBs an d 20 sec for long-duration 
GRBs dZhang fc MeszaroslEiol) . 

Figure[3]plots typical examples of the characteristic am- 
plitudes for neutrino-driven GRBs. The solid line shows the 
amplitude of GWs from the GRB located at the Galac- 
tic centre (r =10kpc), while the short-dashed lines repre- 
sent the GRBs located at r =lMpc (upper) and lOOMpc 
(lower). Here, we assume the total isotropic neutrino energy, 
E = 10 J erg JPopham et al]|l999h . in order to set an upper 
bound of the contributions from the neutrinos to the GW 
amplitudes. Further, extrapolating the ne utrino luminos- 
ity ev olution for the short-duration GRBs iSetiawan et alJ 
120041) to that for the long-duration GRBs, GW amplitude 
from the long GRB was also calculated by simply setting 
tf — 10 sec (solid line, labeled as "Long"). 

Clearly from figure|21 the dominant part of the GWs ap- 
pears at a low-frequency band and they enter the detectable 
band by future observatories. While the high-frequency part 
of GWs rapidly falls down, the characteristic amplitude 
of low-frequency GW asymptotically approaches constant, 
which implies that Fourier component of GWs h(f) is in- 
versely proportional to the frequency, i.e., h cx 1//. This 
behavior typically arises from a sudden change of the neu- 
trino luminosity and can be deduced fr om the so-called zero- 
frequency limit of GW memory (e.g.. [ Braginskv fc Thornei 
119871: lEpsteinllT97l : iBuonanno et alJl2004l 

The results in figure [3] indicate that the planned space 
interferometers, LISA and practical DECIGO can easily de- 
tect the GWs from a GRB at the Galactic centre, though 
the local GRB rate around the Galaxy is extremely small. 
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Figure 3. The characteristic amplitude of GWs from a single 
GRB at the the Galactic centre (solid line). The "Long" duration 
tf is 10 sec and "Short" is 0.1 sec (long-dashed line). On the 
other hand, the short-dashed lines represent the "Long" duration 
GRB at IMpc (upper) and lOOMpc (lower). To compute the GW 
amplitudes, we specifically set the model parameters to £ = n/2, 



tt/6, £ = 10 52 



erg- 



The GRB rate would be increased as the observed volume is 
increasing. For instance, we might observe a few GRB events 
for a decade within lOOMpc. As shown in the figure, how- 
ever, the amplitude of GWs from such a far GRB is quite 
small, which may be comparable or smaller than that of 
GWs generated in the inflationary epoch. In addition, while 
one may detect such GW by an ultimate version of the space 
interferometer whose sensitivity is limited by the quantum 
noise, the detection of GW memory needs a somewhat dif- 
ferent technique from periodic and burst-like GWs, which 
will be an important subject of GW data analysis. 



4 ESTIMATING THE GWB FROM 
NEUTRINO-DRIVEN GRBS 

We are in a position to discuss the contribution of GWs 
from the neutrino-driven GRBs to the background radiation. 
We assume that all GRBs have same energy and opening- 
angle (0 opcn = tt/6), a nd the viewing-a ngle is randomly dis- 
tributed. According to iPhinnevI i200lf) . the sum of the en- 
ergy densities radiated by a large number of independent 
GRBs at each redshift is given by the density parameter 
Agw(/) = Pc 1 {dpGW / dlog f) as: 



£ 2c 



2tt 2 



dz 



^grb(^) 



Gp c Jo 1 + z 

Xr 2 f{\h + (f)\ 2 ),J 



(l + z)H(z) 



.t f /(l + z) ■ 
-/(!+*) 



(8) 



with -Rgrb(^) being redshift distribution of GRBs. H(z) is 
the Hubble parameter. Here the quantity (|ft+(/)| 2 )^,t f rep- 
resents the averaged value of the amplitude |/i+(/)| 2 over the 



Gravitational Waves from Neutrino-Driven GRBs 5 



viewing-angle £ and the duration tf . Note that the combina- 
tion r 2 / 2 (|/i + (/)| 2 )^ itf becomes independent of the distance 
r, and the explicit dependence on the redshift z is eliminated 
after replacing (f,tf) with (/(l + z),tf/(l + z)). 

In order to quantify the amplitude of GWB, we need 
the redshift distribution of GRBs. For a crude estimate 
of the amplitude of Ogw, we use the redsh ift evolution 
inferred from the _Ep-Iumi nosity correlation (lAmati et alJ 
l2002l:lYonetoku et al.ll2004) : 



-Rgrb(z) oc 



(l + zf- 0±1A 
(l + z )0.4±o.a 



z < 1 
z > 1 



(9) 



As for the normalization factor -Rgrb(O), the local observa- 
tion by BATSE indicates 10~ 7 per year per a galaxy, but, 
the jets are highly aligned and the true GRB rate must be 
larger. Then, the local GRB r ate becomes ~ 10 per year 
per a galaxy iFrail et al.ll200l|) . or equivalently, 



-Rgrb(O) = 25 ~ 250 Gpcyr" 



(10) 



Assuming that the orientation of jets is random, the 
ensemble average (\h+(f)\ 2 )^,tf in equation @ is evaluated 
separately, namely, <|M/)| 2 ) e , if oc (|*| 2 > C {|f (f)| 2 ) tf . For 
the average over the viewing-angle £, one finds 



27rsin£dC- 



(11) 



The integration can be done analytically, but the result 
is rather complicated. For a typical set of parameters 
(#o P en, 0) = (tt/6, 1), this gives 0.154. On the other hand, 
for the average of the GW amplitude over the duration, the 
distribution function for the duration of GRBs, T(tt) must 
be taken into account. Here, we use the data set of dura- 
tions (T90) taken from the BATSE 4B catalog 3 adopting 
a logarithmic binning with Alog 10 if = 0.1. The resultant 
distribution becomes a bimodal distribution over the range 
0.015 sec ^ tf ^ 670 sec. Since the function F(/) is the only 
quantity that explicitly depends on the duration tf, the en- 
semble average of the GW amplitude may be independently 
evaluated as: 



(|r(i)| 2 ) tf = 



J\T(t)fT(tf)dt c 



(12) 



jT(t f )dtf 

Figure |1] shows the density parameter of GWB gener- 
ated from neutrino-driven GRBs (solid line), together with 
the sensitivity curves for future missions. To plot the result, 
we adopt the flat cosmological model with the density pa- 
rameters Qm — 0.3 and Q,a = 0.7. The Hubble parameter is 
given by 



H(z) = JW^Af(l + .z) 3 +nA 



(13) 



with the present Hubble parameter Ho being 72 km 
s _1 /Mpc. Note that when evaluating the expression JSJ, the 
range of the integral is restricted to ^ z ^ 20, since most 
of the contribution (~ 80%) to the integral comes from the 
low redshift GRBs z < 4. 

Figure 0] shows that the GWB is broadly distributed 
over the low-frequency band, 10 Hz < / < 10 1 Hz with 
amplitude /iq^gw ~ 10~ 2(} and has almost the flat spec- 
trum with small bimodal bump. To see the contribution from 




10 10" 

f(Hz) 

Figure 4. The energy density parameter of GWB. The solid 
line represents the contribution from neutrinos of all GRBs. The 
dotted line represents one from "Long" -duration GRBs (tf > 2 
sec), and the long-dashed line one from "Short" -duration GRBs 
(if < 2 sec). The dotted-dashed line represents the case of 
i?GRB(0) = 25 Gpc- 3 yr- 1 , J R G RB oc (1 + z) 4 - 6 , (1 + z) - 2 ("low 
rate"). On the other hand, the short-dashed line represents the 
case of highest GRB rates, J?grb(0) = 250 Gpc _3 yr~ 1 , Rgrb °c 
(1 + z) 7,4 , (1 + z) 6 in the cases of z < 1 and z > 1, respectively. 
The horizontal dotted line shows the GWB from the inflationary 
epoch constrained by COBE observations. 



short-duration GRBs and long-duration GRBs separately, 
the ensemble average over the duration (1121 is artificially di- 
vided into the short duration (tf ^ 2 sec) and the long dura- 
tion (tf > 2 sec) , the results of which are respectively shown 
in dashed and dotted lines. One then immediately deduces 
that bimodal shape of the GWB spectrum has originated 
from a distribution of duration. This is mainly because the 
low-frequency part of the gravitational waveform shown in 
figure becomes featureless due to the zero- frequency limit 
of the memory effect. Although the present result was de- 
rived from a specific choice of luminosity evolution (J7J, the 
characteristic behavior seen in figure |1] would remain un- 
changed as long as the collapsar model is applicable to both 
short- and long-duration GRBs. 

Figure 2] indicates that the GWB from neutrino-driven 
GRBs is much smaller than the upper limit of GWB gen- 
erated in the inflationary epoch constrained by the COBE 
observations (the horizontal thick 

dotted-line) and also below the detection limit of fu- 
ture space missions. Note that this tiny amplitude can be 
also deduced from an extrapolation of the discussion by 
iBuonanno et alJ j2004T l based on the following points; i) the 
GRB rate is roughly four orders of magnitude smaller than 
the rate of supernovae and ii) the total released energy by 
neutrinos used in t his paper is thirty tim es smaller than 
theirs. According to IBuonanno et all (I2004T) . the maximum 
amplitude of GWB from cosmological supernovae is esti- 
mated as ~ 10 -12 (see Fig. 4 of their paper). Since 



1 |http://www. batse.msfc. nasa.gov/batse/grb/ca talog/4b /4br_durat ionfflt<mt ensity parameter f^GW is proportional to Rgrb(0)E 2 , 
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(see Eqs.@ and ©), we roughly obtain w 10~ 19 in 

the cosmological GRB case. In this sense, small amplitude 
of GWB may be a natural outcome from the rare even rate 
of GRB. On the other hand, nearly flat spectrum of GWB in 
figure 2] cannot be simply explained by an argument based 
on the the supernovae case, which may be a unique charac- 
teristic of the GRB case. 

As examined in figure 0] (thick-dotted and dot-dashed 
lines), the uncertainty of the normalization in the GRB rate 
11011 does not change the conclusion. In other words, the 
GWB generated from a cosmological population of GRBs 
does not become a preventer of the primordial GWB gener- 
ated in the inflationary epoch. 



5 SUMMARY AND DISCUSSION 

We have discussed the GWs from axisymmetric relativistic 
jets of neutrino-driven GRBs. In order to estimate the am- 
plitude of the GWs, we derived the analytic formulae for 
the gravitational waveform from (ultra-)relativistic energy 
flow. The analytic expression given in equation @ with 
and © depends on the viewing-angle £, the opening-angle 
of the jet # pen, and the velocity of the flow j3, as well as 
the time evolution of luminosity F(t) (or T(f)). Following 
the fitting result obtained from the state-of-the-art numer- 
ical simulations of collapsars, the amplitudes of the GWs 
from a single GRB are estimated and are found that within 
the detection limits for the space-based laser interferome- 
ters like LISA and DECIGO/BBO, if a GRB occurs at the 
region within IMpc. Since the released energy in the jet re- 
gion by neutrinos can be greater than that by matter, the 
low-frequency GWs generated by neutrinos seem to domi- 
nate the one by matter studied bv lSago et al.l l|2004Fl . Note 
that, even if the amplitudes of the GW memory overcome 
the sensitivity curves of interferometers in the spectra, this 
does not directly imply the detection of the signals. The 
detectability should be discussed with a specific data anal- 
ysis technique for the memory effect. On the other hand, 
the GWB from a cosmological population of GRBs are ob- 
tained by summing up the individual GW of the GRBs, 
which turns out to be sufficiently smaller than those from 
the inflati onary epoch and a cos mological population of su- 
pernovae feuonanno et al 

While the sufficiently small amplitude of the GWB 
would be true, there remain several uncertainties in predict- 
ing a precise waveform of the GW from a single GRB jet. In 
particular, when estimating the GW from a long GRB, we 
have extrapolated to use the luminosity evolution of neu- 
trinos for a short GRB. Since the efficient mechanism to 
trigger the long burst is still under debate, one should con- 
tinue to check the present GRB model in more quantitative 
way. As a prelude to more realistic luminosity evolutions 
in the collapsar models, which requires general relativistic 
multidimensional radiation hydrodynamic simulations, our 
findings obtained in a semi-analytic manner should be the 
very first step towards the predictions of the GWs from the 
collapsars and the following formation of GRBs. 



ACKNOWLEDGEMENTS 

We thank S. Yamada and K Sato for informative discus- 
sions. This work was supported in part by the Japan Society 
for Promotion of Science(JSPS) Research Fellowships (TH, 
KK, HK) and a Grant-in-Aid for Scientific Research from 
the JSPS(AT, No.14740157). 



REFERENCES 

Abbott B. et al., 2005, Phys.Rev.D 72, 042002 

Aloy M.A., Janka H .-T., Muller, E., 2004, 

preprint) astro-ph/0408291 ) 
Amati L. et al., 2002, A&A, 390, 81 
Asano K., & Fukuyama T., 2000, ApJ, 531, 949 
Braginsky V.B. & Thorne K.S., 1987, Nature, 327, 123 
Buonanno A., Sigl G., Raffel t G.G., Janka H.-T. & Muller 

E., 2004, |astro-ph/0412000| 
Burrows A., & Hayes J., 1996, Phys.Rev.Lett. 76, 352 
Epstein R., 1978, ApJ, 223, 1037 

Farmer A., J. & Phinney E.S., 2003, MNRAS, 346, 1197 
Fenimore E.E. & Ramirez- Ruiz E., 2000, 

preprint) astro-ph/0004176 ) 
Ferrari V., Matarrese S. & Schneider R., 1999, MNRAS, 

303, 247 

Frail DA. et al., 2001, ApJ, 562, L55 

Fryer C.L., Holz D.E., & Hughes S.A, 2004, ApJ, 609, 288 
Kotake K., Yamada S. & Sato K., submitted to Phys.Rev.D 

Paczynski B., 1990, ApJ, 36 3, 218 

Phinney E.S., 2001, preprint patro-ph/0108028l l 

Proga D., MacFadyen A.I., Armitage P.J., & Begelman 

M.C., 2003, ApJ, 599, L5 
Ruffert M., & Janka H.-T., 1999, A&A, 344, 573 
MacFadyen A.I., & Woosley S.E., 1999, ApJ, 524, 262 
Maggiore M., 2000, Phys. Rep. 331, 283 
Meszaros P., & Rees M.J., 1992, MNRAS, 257, 29P 
Mizuno Y., Yamada S., Koide S., & Shibata K., 2004, ApJ, 

606, 395 

Muller E., Rampp M., Buras R., Janka H.-T., & Shoemaker 

D.H., 2004, ApJ, 603, 221 
Muller E. & Janka H.-T., 1997, A&A, 317, 140 
New K.C.B., 2003, Living Rev. Relativ. 6, 2 
Piran T., 2002, preprint( |g7-qc/0205045| 
Popham R., Woosley S.E., & Fryer C, 1999, ApJ, 518, 356 
Sago N., Ioka K., Nakamura T. & Yamazaki R., 2004, 

Phys.Rev.D 70, 104012 
Schneider R., Ferrari V., Matarrese S. & Portegies Zwart 

S.F., 2001, MNRAS, 324, 797 
Segalis E.B., & Ori A., 2001, Phys.Rev.D 64, 064018 
Setiawan S., Ruffelt M., Janka H.-Th., 2004, MNRAS, 352, 

753 

Seto N., Kawamura S. & Nakamura T., 2001, 

Phys.Rev.Lett. 87, 221103 
Takahashi R. & Nakamura T., 2005, Prog.Theor.Phys., 

113, 64 

Takiwaki T., Kotake K, Nagataki S., & Sato K, 2004, 

ApJ, 616, 1086 
Yonetoku D., Murakami T., Nakamura T., Yamazaki R., 

Inoue A.K. & Ioka K, 2004, ApJ, 609, 935 
Zhang B. & Meszaros P., 2004, Int. J. Mod. Phys. 19, 2385 



